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We report investigations of conductance fluctuations (noise) in doped silicon at low temperatures 
(T< 20K) as it is tuned through the metal-insulator transition (MIT). The scaled magnitude of noise, 
jh, increases with decrease in T following an approximate power law jh ~ T^ 13 . At low T, 7h 
diverges as n/n c crosses 1 from the metallic side. We find that the distribution function and second 
spectrum of the fluctuations show strong non-gaussian behavior below 20K as n/n c decreases through 
1. In particular, the observed distribution function which is gaussian for n/n c >> 1, develops a 
log-normal tail as the transition is approached from the metallic side and eventually it dominates 
in the critical region. 

The Anderson-Mott metal-insulator transition (MIT) has been one of the most studied areas in condensed mat- 
ter (l]J^]. However, even after several decades of work O, there remain many unanswered questions. A fundamental 
issue of considerable theoretical interest is that of conductance fluctuations. The two main aspects of this issue are 
the divergence of noise near the critical region and the onset of non-gaussian fluctuations. The origin of conductance 
fluctuation in the weak localization regime is dominated mainly by the Universal Conductance fluctuations (UCF), 
and is distinctly different from that in the strongly localized insulating regime (<tt— >o = 0) which show activated 
conductance. There are some experimental investigations of the magnitude of conductance fluctuations in thin films 
of In x O y |]l8f at room temperature and in Si inversion layers at low temperatures which showed a divergence 
near the MIT. However, no experimental investigations have looked into the question of non-gaussian fluctuations at 
the MIT. In this report we present the results of the first experimental investigations of non-gaussian conductance 
fluctuations(/ < 10Hz) near the MIT at low temperatures. We have investigated the above aspects of conductance 
fluctuations in bulk single crystals of doped Si as it is driven through the Anderson-Mott transition. Specifically, we 
find that the fluctuations, which have a gaussian distribution in the metallic regime becomes strongly non-gaussian 
with a log-normal distribution at the transition region. This also depends strongly on temperature. We enumerate 
below the important results: 

1. In a 3-D system like doped Si ( Si:P(B)) the noise magnitude shows clear signatures of divergence as n/n c 
decreases through 1 and T — ► 0. 

2. As n/n c — > 1 and the system crosses over to the insulating regime, the conductance fluctuations becomes 
strongly non-gaussian as evident both from the second spectra (defined later) and Probability Distribution 
Function (PDF). 

3. The PDF as well as the second spectra also depend on the temperature and both of them becomes more 
non-gaussian as T is decreased. 

The samples (P doped and in some cases B compensated single crystals of Si) were taken from the same batch 
that were used earlier for transport studies near the MIT M with a range of doping levels n from 3.3xl0 18 /cc to 
1 x 10 19 /cc, covering both sides of the critical doping level (n c ) for MIT (see table, figure 1). They were initially 
characterised by conductivity (a) measurements down to T = O.IK. Figure 1 shows the variation of a as a function 
of T 1 / 2 . When fitted to the simple relation a = uq + mT 1 / 2 , samples with <jq > are taken as metallic. In the 
samples with n/n c < 1 the conductivity for T < 10K follows Mott's variable range hopping (VRH) with the relation 
<r(T) = <j c exp{— (Tq/T) 1 / 4 ). Tq — > as n — > n c from the insulating side. 

Conductance fluctuations (noise) was measured as a function of T for 2K < T < 300K using a five-probe ac 
technique [||. The time series of voltage fluctuations were collected directly into a computer using a 16-bit A/D 
card and was digitally processed to obtain the power spectra <SV(/)- A typical power spectra is calculated from a 
time series of 1-5x10 data points at fixed temperatures which was controlled to within ±lmK. Typically the power 
spectra Sy(f) ~1// Q (a ~ 1) and it scales inversely as the sample volume. At low measurement bias Sy(f) scales as 
V 2 . However, beyond certain bias values, deviation from V 2 dependence was seen at lower temperatures, and will be 
discussed later on. 

In all the samples, the noise magnitude increases as the temperature is decreased below 100K. For T > 100K 
the noise increases with T again and this arises from more classical mechanisms |],[l(J which are outside the scope 
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of this paper. Figure 2 shows a plot of -^p- = / df.S v (f)/V 2 integrated over the experimental bandwidth for a 

representative insulating sample K139 (n/n c =0.94) for T<10K. When measured with low enough biasing field, ^ 
shows an approximate inverse power-law dependence ~ A/T 13 with temperature. The behavior presented in figure 2 
is typical for all the samples on both sides of the MIT, right down to the lowest measured T. 

At lower temperatures in all the samples, the measured noise is sensitive to the biasing field used for the measurement 
as can be seen in figure 2. It shows how the T dependence can change drastically when measured with different 
biases. As T decreases, the bias dependence of - becomes more severe. On increasing the measuring bias beyond 

a characteristic electric field E*(T) (see figure 2 inset), - g2 ' decreases drastically as a function of the measuring bias 
following an inverse power-law. E*(T) in all samples decreased with T. We find that sample heating can explain only 
a small part of this noise suppression. The bias dependence noise has a strong impact on the measured temperature 
dependence of noise, and is clearly seen in all samples. Unless measured with a bias E < E* , the noise tends to "flatten 
out" or saturate at lower temperatures, as seen in figure 2. The cause for this effect is currently under investigation 
and will not be elaborated any further in this paper. All the results presented henceforth have been taken with biasing 
fields E < E* such that Sv(f) oc V 2 and has no bias dependence. 

Figure 3 shows a plot of the Hooge parameter jll]] jjj (=f.NSv(f)/V 2 ) as a function of n/n c for different temper- 
atures for / = 1Hz, N being the total number of carriers in the measured volume. For n/n c > 1 (metallic regime), 
jn is comparatively low. As the doping level decreases and n/n c — > 1 from the metallic side, jjj increases slowly. For 
n/n c ks 1, jjj starts increasing rapidly and at the MIT, changes drastically over several orders of magnitude as it goes 
over to the insulating side. On the insulating side, the observed "fn depends strongly on the temperature (see inset 
3a). At higher temperatures (T > 6K), the value of 7# passes through a sharp peak at a value of n/n c very close to 
but less than 1. It then drops sharply in the insulating side | fl2f . Interestingly this is roughly the same temperature 
(T « 10K ) at which the resistivity starts to become activated and follows the VRH law. 

The inset (b) of figure 3 shows the variation of (3 as a function of the doping level in the system. f3 also has a 
strong n dependence close to the MIT. At the highest doping levels, f3 w 0.5, slowly increasing as the concentration 
is decreased and shows a rather sharp increase while it passes through the MIT to reach a value of 6 in the sample 
with the lowest doping level measured. The varitaion of (3 in the metallic regime as the MIT is approached has been 



discussed elsewhere |13 



On the weakly localized metallic side (n/n c > 2), we have shown earlier jl4|,[15| that the predominant contribution 
to conductance fluctuations arise from the mechanism of UCF p6fl . In this region, the measured noise was found to 
be mostly gaussian (both PDF and the second spectra). However, there is strong departure from such a gaussian 
behavior as n is reduced through n c . A build up of non-gaussian fluctuations as n/n c is reduced through 1 at low 
temperaure (T = 4.2K) can be seen from the probability distribution function (PDF). This is obtained from the 
time series and is the probability of occurence P(AV) of the voltage jump of magnitude AV when the conductivity 
fluctuates in a current biased sample. In all the samples (spanning different regions of n/n c ) the PDF calculated from 
the time series of the conductance fluctuations can be expressed as a sum of a gaussian and a log-normal distribution: 

P(AV) = A g e H ^°»> +f^e 7 k — ^ — ; . A g &n&A ln are the relative weights of the two contributions. Figure 4a 
and b show the PDF of 3 representative samples at 4.2K (in a liquid He bath). For the metallic sample (K348) the 
distrubution is completely gaussian with Ai n =0. However, close to the the transition for the sample K242 one can 
see a gaussian distribution with a log-normal tail. For the insulating sample K240, there is no gaussian component 
(A g = 0) and the PDF is a log-normal distribution. The variation in the relative weights of the two contributions, 
as a function of n/n c , can be seen from the integrated area under the two distributions (Ag for the gaussian part 
and A^at for the log-normal part) normalized by the total area under the PDF (A tot ai) in figures 4(d)and (e). In 
figure 4(d) it can be seen that at T = 4.2K, Ao/A^tai — 1 away from the transition but close to the t ransition 
Ac/ Atotai — 0. A complimentary behavior is seen for Aln / A to tai ■ Figure 4(e) (done for sample K240) shows that 
as T is reduced from the high temperature the gaussian component is reduced while the log-normal component takes 
over. 

Another important measure of non-gaussianity is the second spectra, often referred to as the " noise in noise" . If 
the fluctuations are uncorrelated, then the second spectrum || defined as S^ 2 '(f) — {(AS^ (f)) 2 )um, e should have 
a frequency independent "white" frequency spectrum. S^ 2 \f) is a measure of higher order cumulants. A departure 
of S^ 2 \f) from a white spectrum will be an indication of non-gaussian spectra and finite higher order cumulants. 
Figure 5a shows the second spectrum of the samples at T=4.2K. The spectra (/ < 1Hz) fits the relation S^ 2 '(f) oc f~ p . 
As a function of n/n c , p changes (see inset, figure 5a) from a small value of ~ 0.5 for n/n c > 1 (for white spectrum, 
p = 0)to ~ 1.5 as the insulating state is reached, changing rapidly at n = n c . Figure 5b shows how p varies from a small 
value ~ 0.5 to a large valuer 1.3 in the insulating sample K240 as the temperature is lowered. These observations 
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clearly go to show that as the MIT is approached the dynamics that gives rise to the conductance fluctuation becomes 
increasingly non-gaussian and this is more so as the temperature is lowered. 

The divergence of the conductance fluctuation at the MIT has been predicted in some models. Since the conductance 
a variable with large random fluctuations, it has an important impact on the complete PDF. Some models flTj] propose 
that in the large conductance limit (metallic side) , the relative fluctuations arc smaller and hence the higher cumulants 
of the fluctuations are negligible. This leads to a roughly gaussian PDF. With decrease in conductance, the relative 
fluctuations are larger and the higher cumulants start contributing. On the insulating side, this shows up as a 
log-normal PDF. Since the conductance decreases both as n/n c decreases through 1 as well as T decreases on the 
insulating side, our results are in excellent agreement with this model. However, there is one more important aspect 
in the non-gaussian fluctuation. The departure of the second spectra from a frequency independent white spectrum 
at low temperatures in samples with n/n c < 1 clearly indicates that the non-gaussian fluctuation has a correlated 
dynamics associated with it. We have no clear interpretation as to what gives rise to it but we suggest that this may 
be an indication of a glass like freezing similar to what has been seen in spin glasses ||. Recent theories j| have 
suggested that in the MIT region there is a glass like freezing of the charge system. Analysis of conductivity data 
close to the MIT using various scaling forms has shown that a comprehensive scaling description could be obtained if 
a glass-like transition is invoked at the MIT. Recent noise experiments on 2-D electron system in Si inversion layers || 
has shown that close to the MIT there is an appearance of large low-frequency dynamics which makes the noise blow 
up near the transition. This has been interpreted as a signature of freezing into electron glass. The divergence of 
noise at the MIT is unlikely to arise from the usual dynamics of atomic defect relaxations [6| JTl|Jl^| as the strongest 
coupl ing mechanism that can produce conductivity fluctuation from atomic dynamics is UCF which itself saturates 
at a variance of e 2 /h (Mechanisms like UCF have limits to the fluctuations they can produce, i.e. when UCF is 
saturated). The large jh seen in the insulating samples need much stronger sources of fluctuations. As the system 
goes through the MIT the dopant concentration is reduced and there is no reason why the defect dynamics should 
become more correlated at the MIT (for a n ~ 3 x 10 18 cm -3 the average defect separation is w lOnm). In view of this 
it is most likely that the fluctuation that we see close to the MIT at low temperatures is electronic in origin and is a 
manifestation of the slowing down of the time scale of charge fluctuations as the MIT is approached. On approach 
to the MIT the single particle density of states can show both temporal and spatial fluctuations and if the dynamics 
is glassy then the relaxation time becomes an exponential function of the correlation length Q . We propose that the 
large rise in the low frequency noise at the approach to transition is a manifestation of low frequency dynamics of 
the fluctuation in local density of states which are strongly correlated by coulomb interaction in the absence of strong 
screening. 
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Figure Captions 

Figure 1: Conductivity(a) vs. temperature of the samples Si:P(B). a of the 3 most metallic samples PS24,PS41 and D150 have 
been scaled down by a factor of 5 for convenience. The table shows the relevent sample properties. 

Figure 2: - vs. temperature in the insulating sample K139 measured at different biasing voltages. When measured with 
the lowest bias, the noise scales as T~ f3 for all the samples (for K139, /3 = 6). The inset shows how (AG 2 ) deviates from a G 2 
dependence on application of larger measuring biases. 

Figure 3: The scaled noise jh for different values of n/n c . Inset (a) shows the same graph, magnified around n/n c — 1. Inset 
(b) shows the n/n c dependence of f3. 

Figure 4: 4(a), (b) and (c) shows the PDFs at 4.2K for K348, K242 and K240. 4(d) shows the relative contribution of gaussian 
and lognormal parts (see text) of the PDF as a function of n. 4(e) shows the same as a function of T in sample K240. 

Figure 5: 5(a) Second spectrum S 1 ' 2 ^/) at 4.2K (see text). The graphs have been shifted vertically for visual clarity. The 
inset shows the sharp growth of frequency dependence at the critical concentration n = n c at 4.2K. 5(b) T dependence of the 
exponent p of the second spectra in K 240. 
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